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Abstract 

Inclusive D*^ production is measured in deep-inelastic ep scattering at HERA with the HI 
detector. In addition, the production of dijets in events with a D*^ meson is investigated. 
The analysis covers values of photon virtuality 2 < < 100 GeV^ and of inelasticity 
0.05 < y < 0.7. Differential cross sections are measured as a function of and x 
and of various D*^ meson and jet observables. Within the experimental and theoretical 
uncertainties all measured cross sections are found to be adequately described by next- 
to-leading order (NLO) QCD calculations, based on the photon-gluon fusion process and 
DGLAP evolution, without the need for an additional resolved component of the photon 
beyond what is included at NLO. A reasonable description of the data is also achieved by 
a prediction based on the CCFM evolution of partons involving the /cx-unintegrated gluon 
distribution of the proton. 
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1 Introduction 



Charm quark production in deep-inelastic ep collisions at HERA is of particular interest for 
testing calculations in the framework of perturbative quantum chromodynamics (pQCD). This 
process has the special feature of two hard scales: the photon virtuality and the charm quark 
mass. In the case of jet production the transverse energy of the jets provides a further hard 
scale. In leading order (LO) QCD, the photon-gluon fusion process 7(7 cc is the dominant 
production mechanism. 

Results on inclusive D*^ meson production in deep-inelastic scattering (DIS) have been 
published by the HI and ZEUS collaborations [1-7]. The analysis described in this paper uses 
data collected during 1999 and 2000, corresponding to a larger integrated luminosity of 47 pb^^ 
than used in previous HI publications [1,4,5]. As a result, the production of inclusive D*^ 
mesons is measured in DIS with increased precision. 

The mechanism of charm production is further explored by studying the production of dijets 
in events with a D*^ meson for the first time in DIS. This study is referred to in the following 
as the analysis of "_D*^ mesons with dijets". One of the jets typically contains the D*^ meson. 
According to Monte Carlo simulation studies the jet containing the D*^ meson provides a very 
good approximation of the kinematics of the associated charm quark. The other jet usually 
also gives a good approximation of the energy and direction of the second charm quark or of 
a radiated gluon. Regarding the theoretical predictions, the measurement of jets has a reduced 
sensitivity to fragmentation uncertainties compared to the method based solely on measuring 
the D*^ meson. Thus it corresponds more closely to a measurement of the underlying partons 
and is therefore expected to lead to more reliable theoretical predictions. 

The photon-gluon fusion process of charm quark pair production provides sensitivity to the 
gluon distribution in the proton. The dijet event sample is used to measure the observed gluon 
momentum fraction x'^^. The azimuthal angular correlation A</) between the two leading jets is 
investigated because of its sensitivity to initial state gluon emissions. These measurements are 
compared with QCD calculations based on either coUinear or /cx-factorisation and using gluon 
densities obtained from QCD fits to HERA inclusive DIS data. The photon, in addition to its 
coupling as a point-like object in the hard scattering process, exhibits a partonic structure [8] 
which can be resolved by the hard scale present in the process, and which is described by a 
photon structure function. Dijets are used to measure the observed fractional momentum of the 
parton from the photon taking part in the hard interaction. Comparing this distribution with 
pQCD calculations allows a test for resolved contributions of the photon which go beyond what 
is already included in calculations of the photon-gluon fusion process at next-to-leading order 
(NLO). 

In this paper measurements of single and double differential cross sections for the produc- 
tion of D*^ mesons and D*^ mesons with dijets are reported. They are compared to perturbative 
QCD calculations using different implementations of the evolution of the gluon from the proton. 

2 QCD Models 

QCD models for data corrections and for comparison with measured cross sections are in- 
troduced in the following two sections. Relevant parameters used in the Monte Carlo (MC) 



4 



simulations and the NLO calculations are described and summarized in tabled] 

2.1 QCD Models for Data Corrections 

Monte Carlo models are used to generate charm events and to simulate detector effects in order 
to determine the acceptance and the efficiency for selecting events with a D*"^ meson only 
and with dijets and to estimate the systematic uncertainties associated with the measurements. 
All the events are passed through a detailed simulation of the detector response based on the 
GEANT simulation program [9] and are reconstructed using the same reconstruction software 
as used for the data. 

The Monte Carlo programs RAPGAP [10] and HERWIG [11] are used to generate in DIS 
the direct process of photon-gluon fusion to a heavy (charm or beauty) quark anti-quark pair, 
where the photon acts as a point-like object. In addition, they allow the simulation of charm 
production via resolved processes, where the photon fluctuates into partons, one of which in- 
teracts with a parton in the proton and the rest produces the photon remnant. Both programs 
use LO matrix elements with massive (massless) charm quarks for the direct (resolved) pro- 
cesses. Parton showers based on DGLAP evolution are used to model higher order QCD ef- 
fects. The masses of the c and b quarks are set to = 1.5 GeV and = 4.75 GeV. In 
RAPGAP the hadronization of partons is performed using the Lund String model as imple- 
mented in PYTHIA [12]. For the longitudinal fragmentation of the charm quark into the D*"^ 
meson the Bowler parametrisation [13] is taken with parameters as obtained by BELLE [14]. 
The fragmentation fraction /(c D*+) = 0.257 ± 0.015 ± 0.008 [15] is used. RAPGAP is 
interfaced to HERACLES [16] in order to simulate the radiation of a photon from the incoming 
or outgoing lepton including virtual effects. The simulation of such effects is only available 
for direct processes. For resolved processes, similar QED effects are assumed. RAPGAP is 
normally used for the determination of the detector acceptance and efficiency. The effect of a 
different model on the detector acceptance and efficiency is investigated by using the HERWIG 
program, which is based on the cluster hadronization model. For both models small differences 
in the spectrum of the transverse momentum of the D*"^ meson in comparison with the data are 
corrected by reweighting the spectrum to that observed in the data. 

2.2 QCD IModels used for Comparing with Data 

In this paper the experimental results are compared with predictions considering three active 
flavours (u, d, s) in the proton (fixed-flavour-number scheme FENS) and massive charm quarks 
produced via photon-gluon fusion. The results are also compared with a calculation based 
on the zero-mass variable-flavour-number scheme (ZM-VFNS), where the charm quark occurs 
also as an incoming parton and is treated as massless. In the case of the FENS two different 
pQCD approaches are used: an NLO calculation [17-19] based on the conventional coUinear 
factorisation and the DGLAP evolution equations [20] and another prediction based on k^- 
factorisation and parton evolution according to the CCFM equations [21]. A beauty contribution 
of 1.5 ± 0.5 times the QCD prediction is added to the charm expectations of the calculations 
with massive charm quarks to encompass the range of beauty cross section measurements [22]. 



5 



According to Monte Carlo studies, the contribution of beauty quarks is approximately 3% for 
events with D*^ mesons and 7% for events with D*^ mesons with dijets. The basic parameter 
choices for the various pQCD programs and the range of their variations are summarised in 
table [TJ Each of the variations is performed independently to determine specific cross section 
uncertainties. These uncertainties and the assumed error on the beauty contribution are added 
in quadrature to obtain the total theoretical uncertainty which is shown in the figures as a band. 
The programs implementing the above approaches are discussed in the following. 
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Table 1 : Parton density functions (PDFs) and parameters used in the IVIonte Carlo simulations 
and the NLO programs. The renormalisation and factorisation scales are set equal, /i = = /if 
(apart from CASCADE where /if is kept fixed at /if = a/7T~Q^, where the invariant mass 
squared and the transverse momentum squared of the cc-pair are denoted by s and Q\, respec- 
tively), rric is the charm quark mass and ec and a are the fragmentation parameters according 
to the Peterson and the Kartvelishvili parametrisations. The range of variation for the different 
parameters is also indicated. Each of the variations is performed independently. 

The NLO C(a^) QCD FENS predictions are calculated using the program HVQDIS. The 
CTEQ5E3 [28] parton densities of the proton are used. Charm quarks are fragmented in the 
7P center-of-mass frame into D*^ mesons using the Kartvelishvili et al. [31] parametrisation 
for the fragmentation function which, for HVQDIS, yields a better description of the HI data 
[26] than the Peterson [32] parametrisation. "Decays" of beauty quarks to D*^ mesons are 
parametrised by adapting the longitudinal as well as the transverse fragmentation distribution 
from RAPGAP using the Peterson model with = 0.0080. 

In order to compare parton level dijets of HVQDIS with hadron level dijets of the data, 
hadronization corrections have to be applied to the NLO calculations. They are estimated by 
using the RAPGAP and HERWIG Ivlonte Carlo models described in the previous section. Dijets 
are reconstructed at the parton level from the generated quarks and gluons after the parton 

' For data corrections the CTEQ5L [28] parton density function is used for the proton. 
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showering step, using the same jet algorithm and selection cuts as at the hadron and detector 
levels. For each kinematic bin the ratio of the hadron to parton level cross section is calculated. 
The average values from the two Monte Carlo models are taken as hadronization corrections to 
the NLO predictions. They vary typically between —5% and —20% and occasionally amount 
to —40% and up to +40%. Their uncertainty is taken to be half the difference between the 
predictions of the two models. This uncertainty is added in quadrature to the other theoretical 
uncertainties to obtain the total error for the HVQDIS prediction of D*"^ mesons with dijets. 

The predictions based on the CCFM evolution equation are calculated using the CAS- 
CADE [33] program. In CASCADE the direct process cc is implemented using off- 
shell matrix elements convoluted with the A;T-unintegrated gluon distribution of the proton. 
The parametrisation set-AO [25] is used for the latter. It has been determined from a fit to F2 
data published by HI [34] and ZEUS [35]. Time-like parton showers off the charm quark and 
anti-quark but not off initial state gluons are implemented. The hadronization of partons is per- 
formed in the same way as described in section 2.1 for RAPGAP. The sensitivity of the cross 
section to the parametrisation used for the longitudinal fragmentation of the charm quark into 
the D*"^ meson is investigated by using the Peterson function with parameters as obtained by 
HERA measurements [26, 27] instead of the Bowler function. Resolved photon processes are 
not implemented in CASCADE. 

An NLO QCD calculation of inclusive hadron production in DIS in the ZM-VFNS has 
recently become available [36]. This calculation treats the charm quark as massless. The con- 
tribution of the partonic subprocesses 75 qg and 7(7 qq for the production of charm and 
beauty quarks is considered at LO. The NLO corrections are large in certain kinematic regions. 
To allow a comparison of the data to the predictions from the "massless" calculation, it is re- 
quired that the D*^ meson has a transverse momentum > 2 GeV in the 7p center-of-mass 
frame. At present, predictions from the "massless" approach exist only for inclusive D*^ meson 
production and not for the production of D*^ mesons with dijets. 



3 HI Detector 

The data presented were collected with the HI detector at HERA in the years 1999 and 2000. 
During this period HERA operated with 27.5 GeV positrons and 920 GeV protons colliding at 
a center-of-mass energy of ^/s = 318 GeV. The data sample used for this analysis amounts to 
an integrated luminosity of £ = 47.0 pb~\ 

A detailed description of the HI detector is given in [37]. Here only the relevant components 
for this analysis are described. The positive 2;-axis of the HI reference frame, which defines the 
forward direction, is given by the proton-beam direction. The scattered positron is identified 
and measured in the SpaCal calorimeter [38], a lead- scintillating fibre calorimeter situated in 
the backward region of the HI detector, covering the polar angular range 153° < 6 < 177.8°. 
The SpaCal also provides time-of-flight information which is used for triggering purposes. Hits 
in the backward drift chamber (BDC) are used to improve the identification of the scattered 
positron and the measurement of its angle. Charged particles emerging from the interaction 
region are measured by the Central Tracking Detector (CTD), which covers a range 20° < 9 < 
160°. The CTD comprises two large cylindrical Central Jet drift Chambers (CJCs) and two 
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^-chambers situated concentrically around the beam-line within a solenoidal magnetic field of 
1.15 T. It also provides triggering information based on track segments measured in the r-0- 
plane of the CJCs, and on the ^-position of the event vertex obtained from the double layers of 
two multi-wire proportional chambers (MWPCs). In the central and forward region the track 
detectors are surrounded by a finely segmented Liquid Argon Calorimeter (LAr). It consists of 
an electromagnetic section with lead absorbers and a hadronic section with steel absorbers and 
covers the range 4° < 9 < 154°. 

The luminosity determination is based on the measurement of the Bethe-Heitler process 
{ep ep7), where the photon is detected in a calorimeter located downstream of the interaction 
point in the positron beam direction. 

4 Event Selection 

At fixed center-of-mass energy, ^/s, the kinematics of the inclusive scattering process ep — > eX 
is determined by any two of the following Lorentz-invariant variables: the Bjorken scaling 
variable x, the inelasticity y, the square of the four-momentum-transfer and the invariant 
mass squared W"^ of the hadronic final state. In this analysis these variables are determined 
from the measurement of the scattered positron energy, E'^, and its polar angle, 9'^, according to 



where Ef, is the incident positron beam energy. 

The analysis covers the kinematic region 2 < < 100 GeV^ and 0.05 < y < 0.7. DIS 
events were triggered by requiring signals from the central drift chambers and the multi-wire 
proportional chambers in coincidence with signals from the scattered positron in the SpaCal. 
The identification and selection of the scattered positron is performed as described in [1]. 

D*"^ mesons are reconstructed using the decay chain D*^ — D^ixf — > A'~7r+7r+ (and c.c), 
where the notation vr^ is used for the slow pion. The three decay tracks are measured in the cen- 
tral track detector. For all tracks particle identification is applied using the measurement of the 
energy loss, dE/dx, in the CJCs. The invariant mass of the K^ti^ system is required to be con- 
sistent with the nominal mass within two standard deviations. The signal is extracted from a 
simultaneous fit to the distribution of Am = mK-K-K — f^K-w of the D*^ meson candidates and of 
the wrong sign combinations {K^'K^)'K'f which provide a good description of the shape of the 
uncorrected background. Further details are described in [1,5,39]. The range of the transverse 
momentum and the pseudorapidity of the D*^ meson is restricted to 1.5 < < 15 GeV and 
\ri\ < 1.5, where and rj are defined in the laboratory frame with 77 = — In tan (|) . From the 
fit a total of 2604 ± 77 D*^ mesons is obtained. 

In order to define the meson with dijets" sample, the fcx-cluster algorithm [40] in its 
inclusive mode is applied to the hadronic final state objects in the Breit frame for events contain- 
ing a D*^ meson candidate. The hadronic objects are built by combining the energy depositions 




(1) 





X = — 
ys 
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in the SpaCal and the LAr calorimeter with the track momenta measured in the tracking sys- 
tem. Such objects have improved energy and angular resolution and are well suited for the low 
transverse momentum jets produced in charm events at HERA. For jets with transverse energy 
3 < E^r^^ < 5 GeV the energy resolution is about 20%. When applying the jet algorithm, 
the four-vector of the reconstructed D*"^ meson is used instead of the four-vectors of its three 
decay particles. The jet algorithm is used with the separation parameter set to unity and using 
the .E-recombination scheme, in which the four-vectors of the hadronic objects are added. 

For the dijet selection, the transverse energies of the leading jets in the Breit frame are 
required to be E^r^^^^"^^ > 4(3) GeV, and their pseudorapidities in the laboratory frame have to 
fulfill —1 < ?7iab^'^ < 2.5. Down to these low jet transverse energies Monte Carlo simulation 
studies show a very good correlation of the parton and reconstructed jet quantities for both jets. 
From a fit to the Am distribution a total of 668 ± 49 -D*^ mesons is obtained for events fulfilling 
the dijet requirements. In about 90% of the events the -D*^ meson belongs to one of the two 
leading jets in agreement with Monte Carlo predictions. 



5 Cross Section Determination and Systematic Errors 

The total visible cross section for inclusive -D*^ meson production in deep-inelastic ep scatter- 
ing, including requirements on the DIS phase space and on the kinematics of the D*^ meson, 
is calculated from the observed number of D*^ mesons, A^^^.i, according to 

a(e+p - e^D^^X) = (2) 

C ■ Br ■ e ■ [1 + drad) 

Here Br refers to the branching ratio Br{D*+ D°7r+)-Br{D° K-'k+) = 0.0258 [41] and 
C to the integrated luminosity. The factor e = 31% corrects for the acceptance loss due to the 
track selection cuts and the detector efficiency and resolution. The QED radiative correction, 
(5rad, amounts to —2%. For differential cross sections the data sample is divided into bins and 
the number of -D*^ mesons is extracted in each bin separately. The visible total and differential 
cross section for the mesons with dijets is defined in a similar way. 

The systematic errors on the cross section measurements are estimated as follows (numbers 
are given for the total visible cross section): 



• The trigger efficiency is monitored using data samples with independent trigger condi- 
tions. Its associated uncertainty is estimated to be 2%. 

• The uncertainty in the track reconstruction efficiency leads to an error of 4.8%. 

• To account for possible imperfections in the description of low transverse momentum 
tracks by the Monte Carlo simulation the requirement for the minimal transverse momen- 
tum of the TTg candidate is varied from 120 MeV to 150 MeV. This leads to a cross section 
change of 4%. 

• An error of 3% due to the uncertainty on the dE/dx measurement used for particle iden- 
tification is estimated. 
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Systematic Uncertainties 



D*^ meson 



D*^ + dijets 



±2% 
± 4.8% 
±3% 
-4% 
± 4.9% 
± 1.8% 
± 1.5% 
± 2.5% 
± 1% 

±4% 
± 1% 

+ 8% + 9% 

- 9% - 11% 



Table 2: Experimental systematic uncertainties on the total visible production cross section for 
inclusive D*"^ mesons and for D*"^ mesons with dijets. 



• The systematic error of the D*"^ signal extraction procedure is estimated to be 4.9%. 

• The estimated uncertainties on the measurement of the scattered positron energy E'^ of 
1% and of its polar angle 9'^ of 1 mrad lead to an error of 1.8%. 

• A model uncertainty of 1% is estimated from the difference in the correction factor for 
the acceptance and efficiency obtained with RAPGAP and HERWIG. 

• An uncertainty of 1.5% is caused by the error in the determination of the luminosity. 

• The uncertainty on the -D*^ and branching ratios contributes an error of 2.5%. 

• For the analysis of D*"^ mesons with dijets additional errors of 4% and 1% are taken into 
account. They arise from the uncertainty in the hadronic energy scale of the LAr (4%) 
and the SpaCal (7%), respectively. 

The systematic uncertainties on the total visible cross section are summarized in table[2l For the 
differential cross section measurements the systematic errors are evaluated separately for each 
bin. All contributions are added in quadrature to obtain the total systematic errors. 

6 Inclusive D*"^ Meson Cross Sections 

The cross section for inclusive D*"^ meson production in the DIS kinematic region 2 < Q'^ < 
100 GeV^ 0.05 <y<0.7 and in the visible D*"^ range 1.5 < Pt < 15 GeV and < 1.5 is 
found to be 

c^vis(e+p ^ e-^D*^X) = 6.99 ± 0.20 (stat.) toU (syst.) nb. (3) 



Trigger efficiency 
Track reconstruction efficiency 
dE/dx measurement 
Description of tt^ tracks with low pT 
D*^ signal extraction 
Measurement errors on E'^ and 9'^ 
Luminosity measurement 
Branching ratio 

Model dependence of acceptance and reconstruction efficiency 
Hadronic energy scale of the LAr Calorimeter 
Hadronic energy scale of the SpaCal Calorimeter 
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A comparison with the predictions from HVQDIS and CASCADE is shown in table [3] The 
models include a small beauty contribution as described in section 2. The predictions of both 
calculations are slightly below the data (~ 12%) but are consistent with the data within errors. 
The quoted theoretical errors include the variations of the scale fi, the charm quark mass nic and 
the fragmentation parameters as indicated in table [B All these contributions to the theoretical 
uncertainty are of roughly similar size. These uncertainties and the assumed error on the beauty 
contribution, which leads to a relatively small cross section error, are all added in quadrature 
to define the total theoretical uncertainty. Use of the MRST2004F3NLO [42] instead of the 
CTEQ5F3 parton densities of the proton for HVQDIS results in a 9% decrease of the cross 
section. If for the fcx-unintegrated gluon distribution of the proton in CASCADE the J2003 
set-1 [43] parametrisation is used, instead of set- AO, the cross section increases by 3%. 







HI data 


HVQDIS 


CASCADE 




> e+D*^X) 


6.99 ± 0.20 (stat.) tl% (syst.) nb 


Q.llt'.f, nb 


6.i9t°:^^ 


o-vis(e+p - 


> e+D*±jj X) 


1.60 ± 0.12 (stat.) to\l (syst.) nb 


1.35tO:l^ nb 


1.65tEJ;JE; nb 


o-vis(e+p - 


^ e+D*^i] X) 


0.228 ± 0.014 (stat.) ± 0.011 (syst.) 


n 991 +0.029 

U.ZZ1_Q Q2Q 


U.ZU ( _o.022 


a'vis(e+p - 


e+D*± X) 



Table 3: Comparison of the cross sections for inclusive D*^ meson production and dijet pro- 
duction in association with a D*^ meson and of their ratio with the predictions from HVQDIS 
and CASCADE. The errors on the predictions include the variations of parameters, as indicated 
in table m and the assumed ±33% error on the beauty contribution added in quadrature. 

In figure [H the single differential cross sections for inclusive D*^ production in the visible 
region are showr@ as a function of the event variables Q"^, x and W and of the D*^ observables 
Pt and rj and of the inelasticity z. The latter is defined asz = P- p/P-q = (E — p^) d* I'^yEe, 
where P, q and p denote the four-momenta of the incoming proton, the exchanged photon and 
the observed D*^ meson, respectively. This quantity is a measure of the fraction of photon 
energy transferred to the D*^ meson in the proton rest frame and it is sensitive to both the 
production mechanism and the c D*^ fragmentation function. The measured cross sections 
shown in figure [D are listed in tables |4] and [5] and are in good agreement with previous mea- 
surements from HI [1]. Double differential cross sections as functions of Q'^ and x are listed in 
table [61 

Figure [H includes the expectations from the HVQDIS and the CASCADE programs. The 
ratio R of the theoretical to the measured cross section is also shown for selected distributions. 
The steep fall of the cross section as a function of and x is described by both HVQDIS 
and CASCADE. There is reasonable agreement between HVQDIS and the data for the different 
single differential cross sections with the exception of the medium values of pt and the region 
?7 > 0, where the measured D*^ meson production cross section is larger than predicted. An 

^The bin averaged cross section is shown at the position of the centre-of-gravity of the cross section in that bin 
as calculated by RAPGAP. 
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excess is also observed at small z, a region correlated with the forward direction (77 > 0). 
Already in [1] indications of an excess observed in the data at large pseudorapidities (0.5 < 77 < 
1.5) and small z with respect to the HVQDIS expectation had been reported. The data presented 
here confirm this excess with better statistical precision and across the whole range of Q^, as 
shown in a detailed study of correlations among the observables in D*^ meson production [39]. 
The predictions from the CASCADE program are found to generally agree better with the data 
than those from HVQDIS. An excess of the data over coUinear NLO predictions at 77 > and 
small z has been also observed in the photoproduction of D*^ mesons [44]. 

In figure [2] the inclusive D*^ production cross section with the additional condition on the 
D*^ meson > 2.0 GeV, in order to be able to compare with the ZM-VFNS predictions, 
is shown. The predictions of CASCADE and particularly the one of the ZM-VFNS approach 
are not able to describe the data at large x and large , while HVQDIS is consistent with the 
data within errors. The expectation of the ZM-VFNS calculation is very close to the one from 
CASCADE for and x, however for the forward direction in 77 it is lower than CASCADE 
and similar to the expectation from HVQDIS. The measured cross sections are listed in table|7] 



7 Production cross sections for D*^ Mesons with Dijets 

The production cross section of D*^ mesons with dijets in the kinematic region 2 < < 
100 GeV^ 0.05 <y < 0.7, in the visible D*^ range 1.5 < Pt < 15 GeV and \^\ < 1.5, and 
with jets having Breit frame transverse energies E^r^^^^'^^ > 4(3) GeV, and laboratory pseudo- 
rapidities —1 < ?7iab^'^ < 2.5, is found to be 

c^vis(e+p ^ e+L>*^jj X) = 1.60 ± 0.12 (stat.) toit (syst.) nb. (4) 

The predictions from HVQDIS and CASCADE are listed in tabled As for the inclusive D*^ 
meson analysis the nominal HVQDIS value is lower (16%) than the data, but there is agreement 
within errors. The prediction by CASCADE agrees well with the data. The uncertainty of the 
two predictions due to fragmentation is much reduced compared to the inclusive D*"^ meson 
case. For CASCADE, also the scale uncertainty is reduced while for HVQDIS it remains the 
dominant error contribution which is found to be even larger than for the inclusive production 
of D*^ mesons. Use of the MRST2004F3NLO instead of the CTEQ5F3 parton densities of 
the proton for HVQDIS results in a 1% increase of the cross section. If for the /cT-unintegrated 
gluon distribution of the proton in CASCADE the J2003 set-1 [43] parametrisation is used, 
instead of set- AO, the cross section increases by 1%. 

The systematic uncertainty is reduced significantly in the ratio of the production cross sec- 
tion for D*"^ mesons with dijets to the one for inclusive -D*^ mesons 

, , .r^JJ = 0.228 ± 0.014 (stat.) ± 0.011 (syst.). (5) 

As shown in table |3] the expectation from the HVQDIS program agrees well with the measured 
ratio while CASCADE predicts a somewhat larger value. 
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In figure [3]the differential cross sections for D*"^ mesons with dijets are presented as func- 
tions of the event variables and x, and of the jet variables E^'"^^ = E}^^^ which is the 
maximum transverse jet energy in the Breit frame and the invariant mass Afjj of the dijet sys- 
tem. The data are compared with the expectations from HVQDIS and CASCADE indicated by 
bands, which indicate the total theoretical uncertainty as described in section 2. The uncertainty 
on rric is responsible for the large error of the HVQDIS prediction for low values of E^^"^. The 
quality of the description is more obvious in the ratio of the predicted to the measured cross 
sections, which is also shown in figure [3l Both HVQDIS and CASCADE describe the steep 
fall of the cross section as and x become large, though CASCADE systematically overesti- 
mates the cross section at high values of and x. The distributions of E^^"^ and Mjj are well 
described by both HVQDIS and CASCADE. 

The absolute difference in azimuthal angle in the Breit frame, A0 = |0jet i — 0jet 2I, is shown 
in figure m as a double differential cross section for two ranges of Q^. In the LO photon-gluon 
fusion process ('yg cc) the two jets are expected to be back-to-back, i.e. A0 = 180°. 
The contributions at A0 < 180° arise primarily from hard gluon emissions and fragmentation 
effects. HVQDIS allows radiation of one hard gluon in NLO, while CASCADE includes the 
radiation of one or more gluons in the parton showering process and it contains transverse 
momentum effects of the exchanged gluons. The measured cross sections for large A0 (bin 
two and three) are well described by both theoretical approaches as expected. In order to be 
sensitive to higher order or k-r effects at smaller A0, and because theoretical uncertainties are 
reduced, the following ratio is defined 

norm ,0 theory / „ J2_data ' 

J dO2dA0 / J dQ'MAfP 

A0(bin2+3) A0(bm2+3) 



This double ratio of theory over data is also shown in figure IH as well as the data points to in- 
dicate the experimental errors. At the lowest values of A0 (bin 1), CASCADE is slightly above 
the data, indicating that the fcx-distribution in the unintegrated gluon density in CASCADE is 
too broad. HVQDIS, on the other hand, underestimates the cross section at the lowest A0, both 
at lower and at higher Q^, indicating that in this approach effects beyond NLO are needed to 
match the data. Similar conclusions were obtained in photoproduction from measurements of 
the azimuthal correlations of a D*"^ meson and a jet not associated to the D*^ meson [44] and 
of dijets in events with a -D*^ meson [45]. The measured cross sections are listed in table[8l 

The jet which contains the -D*^ meson, -D*-jet (DjjE and the other jet (OJ) with the highest 
E^r^^ not containing the D*"^ meson are further investigated. The OJ allows a larger region in 
rapidity to be accessed, the more forward direction, compared to the DJ. In figure [5] the cross 
section is shown as a function of the pseudorapidity of the -D*-jet, the other jet, and of the 
difference in pseudorapidity of the two jets, At] = r/Dj — Voj, all measured in the Breit frame. 
The pseudorapidity distributions are reasonably reproduced by both HVQDIS and CASCADE. 
While the region of small values of | A?7|, which might be expected to be particularly sensitive 
to low-x dynamics [46], is well described, small discrepancies are observed for forward going 

^It is typically one of the two leading jets. However this is not required in which case the only constraint on the 
-D*-jet is due to the D*^ meson kinematic cuts. 
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other jets and for large \Ar]\. These are more clearly seen in the ratio i?nornE which has reduced 
theoretical uncertainty. The measured cross sections are listed in tabled 

To further improve the understanding of the charm production mechanism in DIS, the ob- 
servables x'^^ and x'^^ are investigated. At LO they give the observed fraction of the photon 
momentum carried by the parton involved in the hard subprocess and the observed fraction of 
the proton momentum carried by the gluon, respectively. The determination of both quantities 
involves the partons emerging from the hard subprocess, which are approximated by the -D*-jet 
and the other jet. 



The observable x°^^ is defined as 



obs _ \" -rz/i^-J ' \~ J-Z/wo 

(i?*-p*)had ' 

where E* and are measured in the 7p center-of-mass frame. In the numerator (E* — pi) is 
summed over all particles belonging to the two jets and in the denominator (E* — p*)had is the 
sum over all hadronic final state objects. In figure [6] the single differential cross section for the 
production of dijets with a D*^ meson is shown as a function of x'^^^ and double differentially 
in three bins of Q^. The distribution of x"^*^ peaks close to 1 as expected from direct processes, 
but has significant contributions at lower values. The HVQDIS predictions are in reasonable 
agreement with the measured cross section as a function of and they describe the de- 
pendence of x'^^, indicating that there is no need for an additional resolved photon contribution 
beyond what is already included at NLO. CASCADE also provides a reasonable description. 
The expectation by RAPGAP with direct and resolved contributions is similar to the HVQDIS 
prediction. For > 5 GeV^ the data can be described by multiplying the RAPGAP direct 
contribution by a constant factor, independent of However, for Q"^ < 5 GeV^ the data in- 
dicate that a constant factor would not be sufficient, and only the addition of a resolved photon 
component leads to a good description in LO models based on coUinear factorisation. 

The observable x"^^ is defined as 

s " 2E; • 

The single differential cross section for D*"^ meson and dijet production is displayed in figure|7] 
as a function of x^^ and double differentially in three regions of Q^. The ratio i?norm has re- 
duced theoretical uncertainty and is also shown in figure|7l Both HVQDIS and CASCADE with 
the parameter settings and the parton density functions listed in table [T] describe the depen- 
dence of The sensitivity to recent parton density parametrisations has been investigated 
by comparing with the predictions of HVQDIS using the MRST2004F3NLO parametrisation 
and the parametrisations set-B [25] and J2003 set-1 for the unintegrated gluon density in CAS- 
CADE. The differences for the various PDFs are small compared to the large uncertainties of 
the data. In figure |7] the predictions for i?norm using for example MRST2004F3NLO and J2003 
set-1 are compared to the default expectations using CTEQ5F3 and set- AO with HVQDIS and 
CASCADE, respectively. The measured cross sections in bins of x'^^ and x°^^ are listed in 
tables [lO] and [H 



(E* - plU + {E* - pDoj 



Here all bins are used for the normalisation in contrast to i?nomi- 
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8 Summary 



Measurements of the total and differential cross sections for inclusive D*"^ production in deep- 
inelastic ep scattering are presented. In addition, cross sections are measured for dijets produced 
in events with a meson. The general features of the data are described by the QCD pre- 
dictions in the FFNS as implemented in HVQDIS and CASCADE, which are based on either 
coUinear factorisation and DGLAP evolution or /cx-factorisation and CCFM evolution, respec- 
tively. 

Overall, CASCADE matches the inclusive -D*^ data somewhat better in normalisation and 
shape, and in particular in the positive pseudorapidity region. Furthermore the prediction of a 
calculation in the ZM-VFNS, where the charm quark is treated as massless is confronted with 
the data. It is also found to yield a satisfactory description. 

Cross sections for the production of mesons with dijets in DIS are presented as a 
function of various event and jet kinematic variables. Both HVQDIS and CASCADE give 
reasonable descriptions of the differential cross sections with HVQDIS providing a slightly 
better match to the data. The discrepancy observed between both calculations and data for 
azimuthal differences of the two leading jets A0 below 150° indicates that the /cx-distribution in 
the unintegrated gluon density in CASCADE is too broad, and that in the approach of HVQDIS 
effects beyond NLO are needed to match the data. 

The x"^^ dependence of the cross section is described within the experimental and theoretical 
uncertainties by HVQDIS, indicating that there is no need for an additional resolved photon 
contribution beyond what is already included at NLO. A reasonable description is obtained also 
by CASCADE. In addition, the x'^^ distribution is in agreement with the QCD predictions. This 
confirms that the input gluon distributions to the models, obtained from fits to inclusive data, 
provide a good representation of the charm data. 
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02 [GeV^] davis(e+p e+D*^X)/dQ^ [nbGeV-^] 

[2,4.22] 0.92±0.05l^;2^ 

]4.22, 10] 0.358 ±0.017t|];°26 

]10, 17.8] 0.142±0.009;°;°i2 

]17.8, 31.7] 0.065 ±0.0041^06 

]31.7, 100] 0.0124 ±0.0009l:°:o°}i 



a; do-vis (e+p ^ e+D*^X)/dx [nb] 

[2.8 • 10-^ 0.0002] 11940± 587^^?^^ 

]0.0002, 0.0005] 6311±31lt^^? 

]0.0005, 0.0013] 2176± 106;j^? 

]0.0013, 0.0032] 498 ± 33i^^ 

]0.0032, 0.02] 18.4± 2.3^2:4 



W [GeV] dcrvis(e+p ^ e+D*^X)/dW [nbGeV^i] 

[70,110] 0.043 ±0.0021°:°^^ 

] 110, 150] 0.048 ±0.002;°:°^^ 

] 150, 190] 0.041 ±0.002;°;°°^ 

]190, 230] 0.0291 ±0.0018t[5:S]!]^f 

]230, 270] 0.0153 ±0.0015+°:°°28 



Pt [GeV] dcrvis(e+p ^ e+D*'^X)/dpT [nbGeV^^] 

[1.5,2] 3.3±0.3t°j 

]2, 2.5] 3.16±0.20t^j? 

]2.5, 3.5] 1.92 ±0.081^1? 

]3.5, 5] 0.79 ±0.041°;°^ 

]5, 10] 0.104±0.007t°°J° 



V dcrvis(e+p ^ e+D*^X)/d7] [nb] 

[-1.5,-1] 2.15±0.14l°J5 

]-l, -0.5] 2.43 ±0.131°:^^ 

]-0.5, 0] 2.48±0.14l°j^ 

10, 0.5] 2.56±0.15l°Ji 
]0.5, 1] 2.49±0.15l°J8 

11, 1.5] 1.86±0.16l°:i^ 



Table 4: Differential cross sections for inclusive D*"^ meson production in bins of Q'^, x, W, 
and t]. The first error is statistical and the second is systematic. 
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z dcrvis(e+p ^ e"^D*^X)lAz [nb] 

[0,0.1] 5.9±0.8tJ:8 

10.1, 0.2] 12.2± 0.8ti:3 

10.2, 0.3] 11.9± 0.7j:};^ 

10.3, 0.4] 9.4±0.6j:?:|^ 

]o.4, 0.5] ii.i± o.6j:?:^ 

]0.5, 0.7] 8.0±0.4j:?J 

]0.7, 1] 1.35 ± 0.1 li;5;^i 



Table 5: Differential cross sections for inclusive D*^ meson production in bins of z. The first 
error is statistical and the second is systematic. 



[GeV^] 


X 


dCTvi,s(e+p ^ e^D*^X)lAQ'^Ax [nbGeV^^] 


[2, 4.22] 


[2.51 • 10-^ 5.01 • 10-5] 
]5.01 • 10■^ 0.0001] 
]0.0001, 0.000158] 
]0.000158, 0.000251] 
]0.000251, 0.000501] 


2315 ± 525l^6« 
5432 ± 495^635 
4160 ± 3591-^^^ 
2154 ± 1991^^^ 
741 ± 85^^^ 


]4.22, 10] 


[0.0001, 0.000158] 
]0.000158, 0.000251] 
]0.000251, 0.000501] 
]0.000501, 0.001] 


899±116li3i 
757 ± 80t?6 
413 ± 33t|2 
162 ± 15^}^ 


]10, 17.8] 


[0.000251, 0.000501] 
]0.000501, 0.001] 
]0.001, 0.01] 


206 ± 20t}^ 
77 ± 9t^ 
4.9±0.5t!^:^ 


]17.8, 31.6] 


[0.000251, 0.000501] 
]0.000501, 0.001] 
]0.001, 0.01] 


34 ± 71!^ 
46 ± 4l^ 
3.5±0.4i°:3 


]31.6, 100] 


[0.001, 0.00251] 
]0.00251, 0.01] 


4.1±0.4t°:^ 

o.59±o.o7j:[!:!]^ 



Table 6: Double differential cross sections for inclusive D*^ meson production in bins of 
and X. The first error is statistical and the second is systematic. 
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Q2 [GeV^] davis(e+p ^ e+D*^X)/AQ'' [nbGeV-^] 

[2,4.22] 0.54 ±0.031^05 

]4.22, 10] 0.181 ±0.010;°;°2i 

] 10, 17.8] 0.069 ±0.005i;];J]E]| 

]17.8, 31.7] 0.031 ± 0.003 

]31.7, 100] 0.0058 ±0.0005i;;:ol5Jl^ 



X dcrvis(e+p ^ e+D*^X)/dx [nb] 



[2.8 • 10^^ 0.0002] 
]0.0002, 0.0005] 
]0.0005, 0.0013] 
]0.0013, 0.0032] 
]0.0032, 0.02] 



Py [GeV] dCTvi.s(e+p ^ e+D*^X)/Ap-Y [nbGeV^^] 

[1.5,2] 0.64±0.1ll[;;§i 

]2, 2.5] 1.47±0.12l[!;}^ 

]2.5, 3.5] 1.20 ±0.0610 14 

]3.5, 5] 0.65 ±0.031°;°^ 

]5, 10] 0.093 ±0.0071*5™^ 



V dtTvis(e+p e+D*^X)/dr] [nb] 

[-1.5, -1] 0.92±0.08i;];?^ 

]-l, -0.5] 1.18±0.08in^ 

]-0.5, 0] 1.27 ±0.091^4 

10, 0.5] 1.40±0.09l[!;ie 
]0.5, 1] 1.42±0.09l^Ji 

11, 1.5] 1.03±0.10lo.i5 



Table 7: Differential cross sections for inclusive D*"^ meson production for > 2.0 GeV in 
bins of Q'^, x, px and 7]. The first error is statistical and the second is systematic. 



7916 ± 4071^^° 
3340 ± 1861^^^ 
974 ± 60^}^^ 
205 ± 18^4^ 
7.1±l.lt;j 



21 



0^ [GeV^] dCTvis(e+p ^ e+ D*^iiX)/dQ'^ [nbGeV^^] 





[2, 4.22] 
]4.22, 10] 
]10, 17.8] 
]17.8, 31.7] 
]31.7, 100] 


0.197±0.02ltE;°^7 
0.059 ±0.007tJj°J2 

0.028 ±0.004t°;°g^ 
0.0154 ±0.0022t|^:^|]29 
0.0033 ±0.0005i;5:ooo8 




X 


d(Tvis(e+p ^ e+D*^iiX)/dx [nb] 




[2.8 • 10-5, 0.0002] 
]0.0002, 0.0005] 
]0.0005, 0.0013] 
]0.0013, 0.02] 


2688 ± 2881^5^ 
1255 ± 140t]i4 
438 ± 49t?i 

13.0± 1.9^2:? 




i?™'"' [GeV] 


dfTvis(e+p e+£'*^jjX)/d£'™'' [nbGeV^^] 




[4, 6.5] 
]6.5, 10] 
]10, 20] 


0.32±0.02lg;[]2 
0.138 ±0.014lE];j]?^ 
0.018 ±0.004+[]°J° 




Mjj [GeV] 


dcrvis(e+p ^ e+i:)*±jjX)/dA/jj [nbGeV^^] 




[0, 11] 
]11, 17] 
]17, 25] 
]25, 50] 


0.058 ±o.oo5tjj;;j]| 

0.089 ±0.008j:g;j]}? 
0.026 ±0.005;5];^^4 
0.0023 ±0.0012t°:°°J^ 


[GeV2] 




dcrvis(e+p ^ e+D*'^iiX)/dA(t)dQ^ [nbGeV^^ o-ij 


101 


[57.3, 152] 
]152, 166] 
]166, 180] 


0.020 ±0.003looo2 

0.122 ±0.017ig;|^27 
0.142 ±0.01910°^^^ 


]1G, 100] 


[57.3, 132] 
]132, 166] 
]166, 180] 


0.00079 ±0.000161:^:^^0^^ 
0.0054 ±0.0007t^°°25 

0.0124±0.0018t[5:|]J]]!^ 



Table 8: Differential cross sections for the production of D*"^ mesons with dijets in bins of Q^, 
X, E^^"^, Mjj and double differentially in bins of A0 and Q^. The first error is statistical and the 
second is systematic. 
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VD3 dcrvis(e+p ^ e+D*^iiX)/dr]D3 [nb] 

[-2,0.9] 0.066 ±0.010tEJ;gJ]^ 

]0.9, 1.5] 0.66 ±0.0615^:^7 

]1.5, 2] 0.62±0.07j:°;°^ 

12, 3] 0.45±0.05tj];J^ 

13, 5] 0.058 ±0.014t|];°i^ 



VOJ davis(e+p ^ e+i?*±jjX)/drtoj [nb] 

[-2,0.9] 0.048 iO.OlOtol^g^ 

]0.9, 1.5] 0.43±0.06t|];°^ 

]1.5, 2] 0.65±0.08j;[! |]? 

12, 3] 0.48±0.05tg;J^ 

13, 5] 0.117±0.019t!^ |^2i 



|A?7| davis(e+p ^ e+i?*±jjX)/d| A??! [nb] 

[0,0.3] 0.90±0.12i°;?^ 

10.3, 0.8] 0.82±0.09l° °? 
]0.8, 1.4] 0.64 ± 0.07+°- 

11.4, 2] 0.36±0.05lo.o4 
]2,4] 0.070 ±0.014j;°°24 



Table 9: Differential cross sections for the production of -D*-jet and other jet in bins of ?7dj, 
?7oj and Arj. The first error is statistical and the second is systematic. 



x°^' d(Tvis(e+p ^ e+D*^i}X)/dx°J" [nb] 





[0.2, 0.45] 
]0.45, 0.7] 
]0.7, 1] 


0.59±0.14l°j^ 
0.94 ±0. 161°: J3 
3.5 ±0.21°;^ 


[GeV^] 




do-vis (e+p ^ e+D*^]]X)/Ax°^<'AQ^ [nbGeV^^] 


[2, 5] 


[0.2, 0.7] 
]0.7, 1] 


0.112±0.022l°;^i^^ 

0.35 ± o.o4l^:[J| 


]5, 10] 


[0.2, 0.7] 
]0.7, 1] 


0.015 ±0.009l^;^Jg 
0.146 ±0.019j:;];J]5^ 


]10, 100] 


[0.2, 0.7] 
]0.7, 1] 


0.0034 ±0.0007l°™J^ 
0.0184±0.0017l^;{]Jj^;^ 



Table 10: Differential cross sections for the production of -D*-jet and other jet in bins of x° 
and double differentially in three ranges of Q^. The first error is statistical and the second is 
systematic. 
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g2 [GeV^i 



[2, 5] 



logi 



dcrvis(e+p e+D*^iiX)/dxf" [nb] 



[-3.3, -2.4] 
]-2.4, -2.1] 
]-2.1, -1.8] 
]-1.8, -0.9] 



log 



10 • 



-3.3 
-2.4 
-2.1 
-1.8 



-2.4] 
-2.1] 
-1.8] 
-0.9] 



38 ± 7tt 



108 ± 10 



-10 



61 ± 5tl 
3.6±0.5tJ:? 



do-vis(e+p 



-D*^iiX)/dxf'dQ^ [nbGeV^^] 



5.2±1.4tJ:^ 



12.0± 2.0j:} 



5.8±1.1 



+0.8 



o.36±o.iols;:g^ 



]5, 10] 



-3.3 
-2.4 
-2.1 
-1.8 



-2.4] 
-2.1] 
-1.8] 
-0.9] 



1.14±0.55to:i5 

5.o±i.oti:o 

2.4±0.5t5jj 

-0.016 
-0.064 



0.178 ±0.050 



]1G, 100] 



-3.3 
-2.4 
-2.1 
-1.8 



-2.4] 
-2.1] 
-1.8] 
-0.9] 



0.018 
0.023 
+0.04 
0.04 

0.34 ±0.04^°:°^ 



0.192 ±0.055 
0.49 ± 0.07 



0.018 ±0.004 



+0.005 
-0.004 



Table 1 1 : Differential cross sections for the production of i^*-jet and other jet in bins of 
and double differentially in three ranges of Q^. The first error is statistical and the second is 
systematic. 
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Figure 1: Differential cross sections for inclusive D*"^ meson production as a function of Q^, 
X, W, Pt, r] and z. The inner error bars indicate the statistical errors, and the outer error bars 
show the statistical and systematic uncertainties added in quadrature. The bands for the expec- 
tations of HVQDIS and CASCADE are obtained using the parameter variations as described 
in secion 2. Figures a), b) and c) also present the ratio R = crthcory/o"data for the predictions as 
bands, by taking into account their theoretical uncertainties. The inner error bars of the data 
points at i? = 1 display the relative statistical errors, and the outer error bars show the relative 
statistical and systematic uncertainties added in quadrature. 
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Figure 2: Differential cross sections for inclusive D*^ meson production with, compared to 
figure [B the additional requirement > 2.0 GeV for the D*^ meson in the 7p center-of-mass 
frame as a function of Q'^, x, p-r and rj. The inner error bars indicate the statistical errors, and the 
outer error bars show the statistical and systematic uncertainties added in quadrature. The bands 
for the expectations of ZM-VFNS, a "massless" QCD calculation [30,36] and of HVQDIS and 
CASCADE are obtained using the parameter variations as described in section 2. The ratio R 
is described in the caption of figure [B 
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Figure 3: Differential cross sections for the production of D*^ mesons with dijets as a function 
of , X, i?™^" in the Breit frame and Mjj. The inner error bars indicate the statistical errors, 
and the outer error bars show the statistical and systematic uncertainties added in quadrature. 
The bands for the expectations of HVQDIS and CASCADE are obtained using the parameter 
variations as described in section 2. The ratio R is described in the caption of figure [B 
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Figure 4: Double differential cross sections for the production of D*^ mesons with dijets as 
a function of A0 in the Breit frame for two regions in Q^. The inner error bars indicate the 
statistical errors, and the outer error bars show the statistical and systematic uncertainties added 
in quadrature. The bands for the expectations of HVQDIS and CASCADE are obtained using 
the parameter variations as described in section 2. Also shown is the ratio -Rnorm' for which the 
cross section in the last two bins is used for normalisation (for details see section 7). 
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Figure 5: Differential cross sections for the production of -D*-jet and other jet as a function of 
the pseudorapidity of the -D*-jet and the other jet (OJ) and of the difference in pseudorapidity 
= Vdj — Voj of the two jets. The inner error bars indicate the statistical errors, and the outer 
error bars show the statistical and systematic uncertainties added in quadrature. The bands for 
the expectations of HVQDIS and CASCADE are obtained using the parameter variations as 
described in section 2. In addition, the ratio -Rnorm is shown (for details see section 7). 
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Figure 6: Differential cross sections for the production of -D*-jet and other jet as a function of 
and double differentially in three regions. The inner error bars indicate the statistical errors, 
and the outer error bars show the statistical and systematic uncertainties added in quadrature. 
The bands for the expectations of HVQDIS and CASCADE are obtained using the parameter 
variations as described in section 2. For x'^^ the uncertainty of the hadronization correction is 
the dominant contribution to the total theoretical uncertainty. In the lower plot, the RAPGAP 
prediction, indicating the direct and the sum of direct and resolved contributions, is shown. 
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Figure 7: Differential cross sections for the production of -D*-jet and other jet as a function 
of x'^^ and double differentially in three regions. The inner error bars indicate the statis- 
tical errors, and the outer error bars show the statistical and systematic uncertainties added in 
quadrature. The bands for the expectations of HVQDIS and CASCADE are obtained using the 
parameter variations as described in section 2. The ratio -Rnorm is also shown, separately for 
HVQDIS and CASCADE. For the HVQDIS (CASCADE) band the CTEQ5F3 (set-AO) PDF 
is used, the central predictions being given by the full line. The central values using the PDF 
MRST2004F3NLO (J2003 set-1) are indicated as dotted line. 
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